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ABSTRACT 

We investigate the distribution of particle acceleration sites, independently of the actual acceleration 
mechanism, during plasmoid-dominated, relativistic collisionless magnetic reconnection by analyzing 
the results of a particle-in-cell numerical simulation. The simulation is initiated with Harris-type 
current layers in pair plasma with no guide magnetic field, negligible radiative losses, no initial per¬ 
turbation, and using periodic boundary conditions. We find that the plasmoids develop a robust 
internal structure, with colder dense cores and hotter outer shells, that is recovered after each plas- 
moid merger on a dynamical time scale. We use spacetime diagrams of the reconnection layers to 
probe the evolution of plasmoids, and in this context we investigate the individual particle histories 
for a representative sample of energetic electrons. We distinguish three classes of particle accelera¬ 
tion sites associated with (1) magnetic X-points, (2) regions between merging plasmoids, and (3) the 
trailing edges of accelerating plasmoids. We evaluate the contribution of each class of acceleration 
sites to the final energy distribution of energetic electrons — magnetic X-points dominate at moderate 
energies, and the regions between merging plasmoids dominate at higher energies. We also identify the 
dominant acceleration scenarios, in order of decreasing importance — (1) single acceleration between 
merging plasmoids, (2) single acceleration at a magnetic X-point, and (3) acceleration at a magnetic 
X-point followed by acceleration in a plasmoid. Particle acceleration is absent only in the vicinity of 
stationary plasmoids. The effect of magnetic mirrors due to plasmoid contraction does not appear to 
be significant in relativistic reconnection. 

Subject headings: magnetic reconnection — particle acceleration 


1. INTRODUCTION 

Rapid bursts of non-thermal radiation are common 
phenomena in high-energy astrophysics (Pietka et al. 
2015). Examples include solar flares (Lin et al. 2003), 
gamma-ray bursts (Piran 2004; Kumar & Zhang 2015), 
flares of blazars (Abdo et al. 2010; Nalewajko 2013), giant 
flares of magnetars (Mereghetti 2008), and gamma-ray 
flares of the Crab Nebula (Abdo et al. 2011; Tavani et al. 
2011). Despite widely different contexts, they pose simi¬ 
lar theoretical questions: what is the energy source, how 
is this energy dissipated, how are particles accelerated, 
and what is the radiation mechanism. We can address 
these questions by studying the fundamental properties 
of plausible physical mechanisms, and comparing the re¬ 
sults with observations. 

Many high-energy astrophysical environments are 
strongly magnetized, with a magnetization parameter 
a = B 2 /(4:7rw) > 1, where B is the magnetic field 
strength and w is the relativistic enthalpy. In such en¬ 
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vironments, magnetic reconnection, operating in the rel¬ 
ativistic regime, can be the dominant dissipation mech¬ 
anism (Zweibel & Yamada 2009; Uzdensky 2011; Kagan 
et al. 2015). Relativistic magnetic reconnection is a very 
complex and non-linear process, difficult to describe an¬ 
alytically or semi-analytically (Blackman & Field 1994; 
Lyutikov & Uzdensky 2003; Lyubarsky 2005; Giannios 
et al. 2009; Nalewajko et al. 2011). A lot of recent 
progress has been made through numerical studies, in 
particular by ab-initio particle-in-cell (PIC) plasma sim¬ 
ulations for collisionless plasmas ( Zenitani & Hoshino 
2001; Jaroschek et al. 2004; Zenitani & Hoshino 2005a,b; 
Daughton & Karimabadi 2007; Bessho & Bhattachar- 
jee 2007; Zenitani & Hoshino 2007, 2008; Lyubarsky 
& Liverts 2008; Jaroschek & Hoshino 2009; Sironi & 
Spitkovsky 2011; Liu et al. 2011; Bessho & Bhattachar- 
jee 2012; Cerutti et al. 2012b, 2013; Kagan et al. 2013; 
Cerutti et al. 2014; Sironi & Spitkovsky 2014; Melzani 
et al. 2014; Guo et al. 2014; Werner et al. 2014; Guo 
et al. 2015; Liu et al. 2015). 

Significant recent progress was motivated by the prob¬ 
lem posed by the Crab Nebula flares — how to exceed 
the theoretical limit for the observed energy of the syn¬ 
chrotron radiation, £ max — 160 M eV(E/B), where E is 
the electric field strength and B is the magnetic field 
strength, due to the radiation reaction (Uzdensky et al. 
2011; Cerutti et al. 2012a). Implementation of the radi¬ 
ation reaction force in the Zeltron PIC code allowed for 
the self-consistent calculation of radiative signatures (due 
to the synchrotron process) from relativistic reconnec¬ 
tion (Cerutti et al. 2013, 2014). But even in the limit of 
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negligible effect of the radiation reaction, radiation pro¬ 
duced during relativistic reconnection has very interest¬ 
ing properties: strong energy-dependent anisotropy (‘ki¬ 
netic beaming’) and rapid variability observed in certain 
directions (Cerutti et al. 2012b). 

In a typical numerical experiment on collisionless mag¬ 
netic reconnection, the initial conditions include a uni¬ 
form Harris-type current layer (Harris 1962; see Kirk & 
Skjaeraasen 2003 for detailed description of such layers 
in the relativistic case). We can consider a configuration 
in which the current layer is in the (x, z) plane and the 
reconnecting magnetic field is along the x axis. Long and 
thin current layers are prone to two types of instabilities: 
(1) a tearing mode (e.g., Drake & Lee 1977) producing 
a series of plasmoids (magnetic O-points) and magnetic 
X-points along x, which subsequently undergo highly dy¬ 
namic evolution by a series of mergers; and (2) a drift- 
kink mode (Zhu & Winglee 1996; Pritchett et al. 1996) 
bending the layer along z. Standard 2.5-dimensional 
simulations (with 2D position space and 3D momen¬ 
tum space) in the (x,y) plane capture only the tear¬ 
ing mode, neglecting the effects of the drift-kink mode 
which tends to inhibit particle acceleration ( Zenitani & 
Hoshino 2005a, 2007). Although it is possible to stabilize 
the drift-kink mode by introducing a guide magnetic field 
(a uniform ^-component) (Zenitani & Hoshino 2005b, 
2008), the results of recent very large three-dimensional 
simulations suggest that even in the absence of the guide 
field the negative effect of the drift-kink mode on non- 
thermal particle acceleration is only temporary and par¬ 
ticle acceleration recovers on sufficiently long time scales 
(Kagan et al. 2013; Sironi & Spitkovsky 2014; Guo et al. 
2014, 2015). 

Several mechanisms of particle acceleration at different 
acceleration sites in relativistic magnetic reconnection 
have been identified (e.g., Hoshino & Lyubarsky 2012). 
The most direct mechanism, that is especially efficient 
at magnetic X-points, is due to the reconnection electric 
field E z induced by reconnecting magnetic field B x in¬ 
flowing along y (Zenitani & Hoshino 2001; Larrabee et al. 
2003; Lyubarsky & Liverts 2008). Another mechanism 
operates between two approaching plasmoids shortly be¬ 
fore they merge — their motion along x induces a strong 
electric field E z of opposite sign to the main reconnection 
electric field (hence called the anti-reconnection electric 
field; Oka et al. 2010, see also Jaroschek et al. 2004). 
In addition, a first-order Fermi mechanism may oper¬ 
ate within the plasmoids due to their contraction (Drake 
et al. 2006), and more generally due to the curvature 
drifts (Dahlin et al. 2014; Guo et al. 2014, 2015). Re¬ 
cent results of very large PIC simulations of relativistic 
magnetic reconnection have demonstrated that power- 
law energy distributions of particles arise quite efficiently 
(Sironi & Spitkovsky 2014; Guo et al. 2014; Werner et al. 
2014; Guo et al. 2015). 

Since at least three acceleration mechanisms have been 
identified in previous numerical studies of magnetic re¬ 
connection, it is necessary to evaluate their relative im¬ 
portance. In this work, we attempt to address this is¬ 
sue by investigating a representative sample of energetic 
electrons (i.e., all individually tracked particles that ex¬ 
ceed a certain energy threshold). We propose a simple 
method of identifying the main acceleration episodes for 
each particle, and use the distribution of these episodes 


to make a crude classification of the main acceleration 
sites. The results presented here are based on a single 
2.5-dimensional PIC simulation of relativistic magnetic 
reconnection in an electron-positron plasma, initiated 
from an unperturbed Harris-type current layer, and per¬ 
formed with the Zeltron code. In Section 2, we present 
the simulation setup. In Section 3, we present the de¬ 
tailed results of our simulation, including the analysis of 
individual histories for a representative sample of high- 
energy particles. A discussion follows in Section 4, and 
conclusions in Section 5. 

2. SIMULATION SETUP 

We perform a 2.5-dimensional PIC simulation of rela¬ 
tivistic magnetic reconnection in pair plasma using the 
Zeltron code 8 (Cerutti et al. 2013) in an (x,y) do¬ 
main with periodic boundaries. The domain size is 
L x = L y = 800pcb where po = m e c 2 /(cBq) is the nominal 
gyroradius and Bq is the characteristic value of the re¬ 
connecting magnetic field strength. The coordinate lim¬ 
its are 0 < x < L x and 0 < y < L y . The initial condition 
is a pair of relativistic Harris-type current sheets (Kirk & 
Skjaeraasen 2003) parallel to each other, oriented along 
the x coordinate (with the current flowing along Ez), and 
centered at y = 0.25 L y and y = 0.75 L y . We set the ini¬ 
tial magnetic field as B x (y) = —Hotanh[(p — 0.25 L y )/S 
for 0 < y < 0.5 L y , and B x (y) = L> 0 tanh (y — 0.75 L y )/S 
for 0.5 L y < y < L y , where S is the characteristic half¬ 
thickness of the initial current layers. We do not have 
any kind of initial perturbation of the magnetic field 
(hence B y = 0), nor any guide magnetic field (hence 
B z = 0). The gradient of magnetic field B x along y is 
supported by the plasma pressure (pressure balance) and 
by the electric current (Ampere’s law), both provided by 
a drifting pair plasma component characterized by ini¬ 
tial dimensionless temperature @d = kT^/(m e c 2 ) = 1 
and by dimensionless drift velocity in the Ez direc¬ 
tion jdd = 0.6. The number density of the drifting 
plasma (including both electrons and positrons) is given 
by n d (y) = nd,o cosh -2 [(y—0.25L y )/ 5] for 0 < y < 0.5 L y , 
and rid(y) = rid,o cosh -2 [(y — 0.75 L y )/5\ for 0.5 L y < 
y < L y , where nd,o = IdB^/(&'KQdim e c 2 ) is the peak 
number density and 7d = (1 — /3^) -1 / 2 = 1.25 is the 
drift Lorentz factor. The temperature and the drift ve¬ 
locity of the drifting electrons and positrons determine 
the thickness of the current layer 5 = 2©dPo/(7dAd) = 
(8/3)po (Kirk & Skjaeraasen 2003). We also have a 
background plasma component of dimensionless temper¬ 
ature ©b = kT^/(m e c 2 ) = 1 and uniform number den¬ 
sity rib corresponding to a background magnetization of 
(To = B^ / {Ai:n\ i m e c 2 ) = 16, and hence to a drifting-to- 
background density ratio of nd,o/^b = 7d cr o/(2@d) = 
10. We do not include the effect of radiation reaction 
from the synchrotron process. The simulation grid has 
N x xN y = 2048 2 cells, so that the cell size is Ax = Ay = 
L x /N x ~ 0.4po — 0.15(5. Initially, each cell contains 256 
particles including background and drifting populations 
of both electrons and positrons. The density structure 
for drifting particles is obtained by adjusting particle 
weights; the weights of background particles are uniform. 
The simulation is run for N t = 10 4 time steps, with 

8 http://benoit.cerutti.free.fr/Zeltron/ 
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each time step given by uioAt = 0.9Ax/\/2po — 0.25, 
where = eL> 0 /(ra e c) = c /po is the nominal cyclotron 
frequency. Therefore, the total duration of the simula¬ 
tion is cjo^max = N t uJoAt ~ 2500, which corresponds to 
ctm&x/L x ~ 3.1 light-crossing times of the simulation do¬ 
main. 

3. RESULTS 

Here, we analyze the evolution of just one of the two re¬ 
connection layers present in our simulation, the one cen¬ 
tered at y — 0.75 L y , as the results are generally similar 
for both layers. Each layer develops a different set of plas- 
moids, the details of which are chaotic and unpredictable 
a priori due to the initial random noise in the particle 
distribution. In addition to (x,y) snapshots (cuts in £), 
we will rely heavily on (x,t) spacetime diagrams (cuts 
in y centered on the current sheet). Also, we will focus 
on the electrons only, as the behavior of the positrons is 
highly symmetric to that of the electrons (they co-move 
along x and counter-move along z). In Section 3.1, we 
describe the evolution of plasmoids in the reconnection 
layer using spacetime diagrams. Then, in Section 3.2, we 
will characterize the acceleration of energetic particles in 
the context of the plasmoid evolution. 

3.1. Evolution of plasmoids 

In Figure 1, we present snapshots of the (x,y) maps 
of electron number density and magnetic field lines for 
several simulation times. Each snapshot shows only a 
fragment of the simulation domain for 0 < x < L x and 
550po < y < 660pce centered on the reconnection layer of 
interest (y = 500po)- As has been demonstrated in nu¬ 
merous similar studies, most of the interesting physical 
activity is concentrated in the vicinity of the reconnec¬ 
tion layers. At t = 0, the system consists of a thin current 
layer and magnetic field B x , both uniform along x. By 
ujot = 249 0.31(L x /po)j the tearing-mode instability 

has disrupted the current layer and produced a chain of 
plasmoids, also known as magnetic islands or magnetic 
flux ropes in 3D, characterized by high particle density 
surrounded by closed magnetic field lines (they are as¬ 
sociated with the peaks of the magnetic vector potential 
A z ). They are separated by magnetic X-points charac¬ 
terized by low particle density (the saddle points of A z ). 
At this early stage, the plasmoids are significantly elon¬ 
gated in the x direction. Subsequently, the plasmoids 
are accelerated along ±x, merge in a hierarchical order, 
and also circularize and maintain a roughly round shape 
after each merger. By uoot = 994 ~ 1.24(L x /po) only 
2 plasmoids remain of the original set, although small 
secondary plasmoids emerge occasionally from the ma¬ 
jor X-point (i.e., the X-point that survives to the end 
of the simulation). Because the x = 0 = L x bound¬ 
ary is periodic, we have the freedom of translating the 
simulation domain along x, and we chose a translation 
(fixed for the entire simulation time) that sets the ma¬ 
jor X-point at the boundary. The location of the major 
X-point cannot be known a priori , as it depends on the 
entire evolutionary history of the plasmoids. By choos¬ 
ing such a translation, we can ensure that no plasmoid 
from the original chain ever crosses the boundary. At 
ujot = 1492 ~ 1 . 87 (L x /po) , a secondary vertical current 
layer (oriented along the ^/-direction) becomes visible be¬ 
tween the two large plasmoids in an early stage of their 


merger. At uj$t = 1740 2.1&{L X /po), the final merger 
is almost complete with a newly formed plasmoid escap¬ 
ing along the —y axis. Plasmoids forming due to the 
tearing mode instability of the secondary vertical current 
sheets were previously seen by Daughton & Karimabadi 
(2007) and Sironi & Spitkovsky (2014). The evolution 
of plasmoids is balanced in the sense that at any given 
instance the existing plasmoids are similar to each other 
and roughly equidistant . 9 

Figure 1 suggests that the evolution of the reconnection 
layer can be projected onto the x axis, preserving most of 
the useful information. Therefore, it is convenient to rep¬ 
resent the evolution of plasmoids on an (x, t) spacetime 
diagram (e.g., Samtaney et al. 2009; Fermo et al. 2011). 
To create a spacetime diagram, for a set of 500 regular 
time steps we integrate the parameters of interest within 
a stripe defined by 590po < y < 610po or Ay = 20po? 
marked with white dotted lines on Figure 1 . At early 
stages of the simulation (uot < 800 = L x /po ) the inte¬ 
gration along y includes whole dense plasmoid cores, and 
at late stages of the simulation only the central portions 
of the plasmoids are included. 

The upper left panel of Figure 2 shows the spacetime 
diagram (x, t) of electron number density n e (both drift¬ 
ing and background combined). Many features that were 
identified in Figure 1 — the onset of the tearing mode 
instability, the hierarchical evolution of plasmoids, the 
emergence of secondary plasmoids, and even the forma¬ 
tion of the secondary vertical current layer shortly be¬ 
fore the final merger — can be seen here at very high 
time resolution. Some details of the plasmoid evolu¬ 
tion are not represented on this diagram — the shape 
of the plasmoids (elongated/circular), or the secondary 
plasmoids produced during the final merger and propa¬ 
gating along —y. On the other hand, the spacetime dia¬ 
gram is very useful for describing the kinematics of plas¬ 
moids. Because of the balanced evolution of plasmoids, 
each merger involves plasmoids of similar size and results 
in a cancellation of most r-momentum. Correspondingly, 
the merged plasmoids are very slow, almost stationary. 
Their bulk acceleration is not very rapid, and they only 
reach trans-relativistic velocities (T/3 ~ 1 ) shortly before 
a merger. Only the small secondary plasmoids appear to 
be propagating with strongly relativistic velocities. It is 
also apparent that the plasmoids attract each other with 
a force that is increasing with decreasing distance; this 
can be understood simply as due to Ampere’s force law. 

The upper right panel of Figure 2 shows the space- 
time diagram of the average electron Lorentz factor ( 7 ). 
This parameter combines the intrinsic energy distribu¬ 
tion of electrons due to their random motions with their 
bulk motion. At t = 0 both the background and drifting 
electrons have Maxwellian distribution with kT = m e c 2 , 
which corresponds to ( 7 ) ~ 3.3. In addition, however, 
the drifting electrons are boosted along the z-axis by drift 
velocity /3d = 0.6, so the actual initial average Lorentz 
factor is ( 7 ) ~ 4.3. For comparison, at the end of the 
simulation the average Lorentz factor of the background 

9 The evolution of plasmoids is very sensitive to the introduc¬ 
tion of initial perturbations. We checked for the effect of long- 
wavelength initial perturbation of the magnetic field. Even at very 
modest perturbation amplitude of 0.03£?o, we observed the rapid 
creation of a major plasmoid resulting in a markedly unbalanced 
subsequent evolution. 
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Figure 1 . Maps in x,y coordinates of the total electron density (color scale ; arbitrary units) and magnetic field lines (cyan lines) taken 
at regular intervals of simulation time indicated on the left axis. Dotted white lines mark the region with Ay = 20po, from which the 
spacetime diagrams were extracted. 


electrons is ( 7 ) ~ 7.0, hence the average background elec¬ 
tron gained a Lorentz factor of (A 7 ) ~ 3.6 = 0.22<Jo. As 
the plasmoids grow, they develop a distinct temperature 
structure (also noticed by Melzani et al. 2014), being 
clearly hotter at the edges (( 7 ) ~ 12 ) than in the center 
(( 7 ) 8 ). A closer comparison with the density diagram 

reveals that high-temperature regions (plasmoid shells ) 
are located outside the high-density plasmoid cores. In 
the spacetime diagram, this structure is seen mainly 
for plasmoids with cores larger than Ay, in which case 
the diagram cuts across them. Full (x, y ) maps of ( 7 ) 
(not shown) reveal that the high-temperature shells com¬ 
pletely surround the high-density cores. This structure 
is very persistent even after multiple plasmoid mergers, 


which indicates that the drifting particles that dominate 
the dense cores do not mix with the background parti¬ 
cles energized in the magnetic X-points that dominate 
the hot shells. The persistency of this structure could 
be specific to 2D simulations, where individual particles 
are constrained by conservation of the canonical momen¬ 
tum P z + ( q/c)A z , where P z is the ^-component of the 
particle momentum, and A z is the ^-component of the 
magnetic vector potential, and hence they are confined 
to a given magnetic field line (e.g., Jones et al. 1998). On 
the spacetime diagram of average Lorentz factor we can 
see additional weak structures induced after each merger 
and propagating far into the low-density regions. Also in 
the low-density regions, we observe fan-like structures of 
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Figure 2. Spacetime diagrams of the reconnection layer averaged over the region with Ay = 20 po. Upper left - electron number density 
(arbitrary units). Upper right - average electron Lorentz factor. Middle left — ^-component of the total electric field. Middle right 
- ^-component of the non-ideal electric field. Lower left - particle locations at the beginning of the main acceleration episodes for a 
representative sample of electrons that exceed Lorentz factor of 20. The color indicates the mean velocity component (v z ) /c during the 
main acceleration episode. Lower right - Regions of the spacetime diagram used for classification of the main acceleration sites: magnetic 
X-points (blue), plasmoid mergers (red), and the plasmoids themselves (green). 
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particles being energized by the electric field in the mag¬ 
netic X-points. However, the average particle energies 
in the low-density regions are significantly lower than at 
the plasmoid edges. High particle energies are also found 
in the secondary plasmoids propagating with relativistic 
velocities. They clearly stand out on the ( 7 ) diagram 
despite having low density contrast. 

The middle left panel of Figure 2 shows the space- 
time diagram of the ^-component of the total electric 
field E z . We note that other components of the elec¬ 
tric field - E x , E y - take much smaller values along 
the reconnection midplane. The middle right panel 
shows the ^-component of the non-ideal electric field 
(E z + (ft;) By - (/ 3 y )B x ), where (ft) = (£f =1 ft,ft/^ 
is the mean value of the i-component of velocity of all 
N particles (electrons and positrons) in each cell of in¬ 
terest. Comparing these two panels, we note that the 
total electric field that determines particle acceleration 
is dominated by the ideal component (more specifically 
by (/3 X ) By) induced by the bulk motions of the plasmoids 
along the x-axis. The non-ideal electric field forms three 
types of distinct structures: regions of negative reconnec¬ 
tion field in the vicinity of the magnetic X-points , 10 re¬ 
gions of positive anti-reconnection field between pairs of 
approaching plasmoids about to merge with each other, 
and rapidly propagating regions of positive field following 
each plasmoid merger on either side of the merged plas¬ 
moid. As we will demonstrate in Section 3.2, these are 
some of the main sites of particle acceleration during our 
simulation, although many particles are also accelerated 
in plasmoids by ideal electric fields. 

The lower panels of Figure 2 are discussed in Section 
3.2.4. 

3.2. Particle acceleration 
3.2.1. Momentum distribution of all particles 

In Figure 3 we show the time evolution of the spec¬ 
trum N(u) of background electrons, where u = 7/3 is the 
dimensionless particle momentum (equivalent to particle 
energy for u 1). The spectrum is compensated by u 2 
in order to reveal at which value of u the particles con¬ 
tain most momentum. The spectrum of the background 
positrons is identical to that of the background electrons. 
The initial momentum distribution of background elec¬ 
trons is the Maxwell-Jiittner distribution for dimension¬ 
less temperature ©b = kT]j/(m e c 2 ) = 1. In time, the 
distribution extends into higher energies, developing an 
apparent power-law section of increasing hardness. Mea¬ 
suring the power-law index a defined by N(u) oc u~ a is 
difficult because the power-law section is relatively short, 
and both the low-momentum and high-momentum cut¬ 
offs are relatively broad (for more discussion see Werner 
et al. 2014). By considering a series of differently com¬ 
pensated spectra (not shown), we find that the final spec¬ 
trum of background electrons (uot = 2486; solid blue 
line in Figure 3) is consistent with a ~ 1.6, which is in 
agreement with the scaling results of Werner et al. (2014) 
(see their Figure 3, noting that here we have <to = 16 

10 The sign of the reconnection electric field E z in magnetic 
X-points is determined by the sign of the reconnecting magnetic 
field gradient dB x /dy, which is opposite for the two layers in the 
simulation. 



Figure 3. Momentum distribution of all background ( solid lines) 
and drifting ( dashed lines ) electrons for 6 simulation times: t — 0 
(black), wot = 497 (brown), wot = 994 (red), wot = 1492 (orange), 
wot = 1989 (green), and wot = 2486 (blue). The energy limit 
7 = 10 for the selection of a representative sample of energetic 
background electrons is indicated with the vertical magenta line. 

and L/aoPo = 50), and also with the results of Sironi & 
Spitkovsky (2014) and Guo et al. (2014, 2015). 

3.2.2. Main acceleration episode 

Our main goal is to investigate the mechanisms of par¬ 
ticle acceleration during the simulation by analyzing a 
representative sample of individual particles with their 
complete histories. The total number of particles in our 
simulation is about 10 9 . It is unfeasible to record detailed 
history of each particle for analysis. Instead, we ran¬ 
domly selected a significant sample of 10 5 tracked back¬ 
ground electrons , for which we recorded their detailed 
history for 500 time steps evenly distributed over the sim¬ 
ulation time. Of this sample, we are interested only in 
the highest-energy particles, and we arbitrarily choose a 
threshold minimum value of the Lorentz factor 7 t hr = 10. 
As is illustrated in Figure 3, this threshold is well above 
the peak of the initial momentum distribution of back¬ 
ground electrons. We selected a representative subsam¬ 
ple of 14169 energetic particles that exceeded 7 thr at least 
once during the simulation, and also at any time crossed 
the reconnection layer under consideration (y = 0 , where 
y = y — 0.75L y ). Using such a representative sample 
of energetic particles is essential to obtain a complete 
picture of acceleration mechanisms at work during the 
simulation. 

For each energetic particle we identify what we call 
the main acceleration episodes. They are defined as the 
shortest mutually exclusive intervals of time [t\ : £ 2 ] dur¬ 
ing which the particle gains energy 7 ^ 2 ) — 7i(U) > 
^( 51 ,max 7 i,min )5 where 7i,max} 7i,min 3,re the global ex¬ 
trema of the particle energy history 7 $(£), and k is a 
constant. Here we adopt k = 0.5 after experiment¬ 
ing with smaller values; hence during the main accel¬ 
eration episode a particle gains energy that is at least 
one half of its total energy range during the simulation. 
A particle may have m > 1 non-overlapping main ac- 
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celeration episodes since 7 (t) is usually non-monotonic. 
In such cases, each main acceleration episode will be 
treated as contributing to the particle maximum energy 
7max with 1/m weight. Each main acceleration episode 
is characterized by the initial location of the particle 
along the reconnection layer x(ti), the episode duration 
At = t 2 — ti, and by the drift displacement along the z 
axis Az = z(^) — z(ti), or equivalently by the mean z 
velocity (v z ) = Az/At. 

3.2.3. Individual particle histories 

Before attempting to characterize this entire subsam¬ 
ple, we will first consider selected individual examples. 
We can roughly classify the energetic particles into 4 
groups, according to the location of their main acceler¬ 
ation site: (1) magnetic X-points, (2) regions between 
merging plasmoids, (3) plasmoids, (4) other/multiple 
sites. Each group is illustrated with 4 individual exam¬ 
ples in Figure 4, and one such example for each group is 
analyzed in detail in Figure 5. 

Particles belonging to group (1) gain most of their en¬ 
ergy in the vicinity of magnetic X-points. A detailed 
analysis of a representative example (‘particle 128’) is 
shown in the first panel of Figure 5. The particle gains 
most energy from 7 ~ 5 at wot = 330 to 7 ~ 39 at 
wot = 700, starting at the major X-point at x ~ 0 and 
drifting towards x ~ 100po- This corresponds well with 
the only main acceleration episode wot G [380 : 580]. At 
the same time, the particle is trapped in the reconnection 
layer (y G [—10 : 10] po) and it drifts monotonically in the 
+z direction on a Speiser-type orbit (Speiser 1965; note 
the oscillations in y and u y ). The particle feels roughly 
constant electric field E z — 0.12H 0 and oscillating mag¬ 
netic field B x peaking at ~ ±0.3F>o (in phase with y). 
One can roughly estimate that such electric field predicts 
a maximum energy gain of A 7 ~ \E Z / Bo\woAt ~ 44, 
which is more than the actual energy gain A 7 = 34 
of ‘particle 128’. The particle z -momentum increases 
until (Jot ~ 600, after which the particle starts inter¬ 
acting with the magnetic field B y (peaking at ~ 1.4F>o 
at wot ~ 750) of the plasmoid centered at x 140po, 
which deflects it towards -\-x. At wot ~ 700, the parti¬ 
cle stops drifting towards -\-z and the acceleration stops. 
During wot G [330 : 700], the particle was displaced 
by Az ~ 310po, which means the average drift velocity 
(p z ) = A z /(wo At) ~ 0.84. Subsequently, the particle 
follows the plasmoid on a wide slightly irregular orbit, 
and it gains additional energy, most likely through cur¬ 
vature drifts (cf. Figure 6 in Oka et al. 2010). 

Group (2) consists of particles that gain most of their 
energy between merging plasmoids. A detailed analy¬ 
sis for a representative example (‘particle 40’) is pre¬ 
sented in the second panel of Figure 5. The parti¬ 
cle gains most energy from 7 ~ 3 at wot = 1530 to 
7 ~ 37 at wot = 1700 between two merging plasmoids 
at x ~ 400po- This is captured by a main accelera¬ 
tion episode at wot G [1555 : 1640]. The particle passes 
through the reconnection layer midplane (y ~ 0) exactly 
at wot = 1530, at which moment it starts drifting sharply 
towards the — z axis. The total displacement of the par¬ 
ticle along z is Az = — 145po, which means the aver¬ 
age drift velocity of (/3 Z ) ~ —0.85. The particle feels 
an anti-reconnection electric field varying in the range 
E z ^ (0.2 — 0.28)F>o• Such electric held predicts a max¬ 


imum energy gain of A 7 ~ 41, slightly higher than the 
actual value of A 7 = 34. Consequently, the particle ends 
up located below the major plasmoid, oscillating in y and 

z. 

Group (3) consists of particles that gain most of their 
energy within plasmoids. A detailed analysis for a rep¬ 
resentative example (‘particle 95’) is shown in the third 
panel of Figure 5. The particle gains energy gradually 
from 7 7 at c o^t = 210 to 7 ~ 50 at wot = 1520, for all 

this time being confined to plasmoids and following them 
in their evolution. Two main acceleration episodes were 
identified at wot G [460 : 770] and wot G [1410 : 1750]. 
As the particle oscillates across the plasmoid, its energy 
fluctuates due to the superposition of the bulk motion 
of the plasmoid. However, during some of the sharp 
turns in u x the particle gains more energy than it lost 
during the previous turn. This is seen most clearly at 
wot ~ 570 and wot ~ 1440. These two turns happen 
at the trailing edge of the plasmoid, where the parti¬ 
cle experiences a strong electric field E z ~ —0.3 B 0 that 
is not balanced by the motional field since u x ~ 0 and 
B x ~ 0. We tentatively identify this acceleration mecha¬ 
nism as the curvature drift (e.g., Dahlin et al. 2014; Guo 
et al. 2015; Li et al. 2015). In contrast to the particles 
accelerated in magnetic X-points (group 1) or between 
merging plasmoids (group 2 ), this particle does not ex¬ 
perience a smooth drift along the z axis. While there 
is a systematic drift towards +z (with the average drift 
velocity (/3 Z ) ~ 0.36), the z-momentum (u z ) oscillates 
with a period twice shorter than oscillations of u x and 
u y . Shortly after the turns in u X: most of the momen¬ 
tum gain is seen in the x coordinate. At the peaks of 
the r-momentum, the particle tends to cancel both its y- 
and z-momentum components. The particle does briefly 
feel the anti-reconnection electric field E z ^ 0.3F>o at 
wot ~ 1620, but the resulting additional acceleration is 
insignificant because of the very short time spent there. 

Group (4) consists of particles for which the accelera¬ 
tion mechanism is different from the 3 discussed so far, 
or there is evidence for more than one acceleration mech¬ 
anism during the particle history. The particular cases 
shown in Fig. 4 include acceleration in the vicinity of 
plasmoids created in a recent merger (particle 113), ac¬ 
celeration between merging plasmoids for a particle pre¬ 
viously bound in one of them (particle 145), and accel¬ 
eration of unbound particles in the trailing edge of a 
plasmoid (particles 263 and 279), showing also that a 
particle accelerated between merging plasmoids can be 
accelerated further by a different mechanism. 

In the fourth panel of Figure 5 we present a de¬ 
tailed analysis of ‘particle 113’, accelerated after a nearby 
plasmoid merger but outside the merger product. The 
main acceleration episode is at wot G [880 : 1080]. At 
wot ~ 780, the particle enters a magnetic X-point at 
x ~ 400po with 7 ~ 5 and undergoes acceleration by the 
reconnection electric field E z ~ —0.1 Ho, being displaced 
by Az ~ 150po- At wot zz 1030 with 7 21 , this mode 

of acceleration stops as the particle approaches a plas¬ 
moid newly formed in a recent merger. The electric field 
switches sign and increases very sharply to E z ~ 0.37Hq, 
together with a strong increase in the magnetic field 
By ~ 1.35Ho« This forces the particle to bounce off the 
plasmoid and accelerate rapidly to 7 ~ 36 at wot = 1100. 
The sharp increase in the electric field is different from 
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particle 128 





particle 113 




particle 210 


particle 133 











particle 273 




Figure 4. Spacetime diagrams of the trajectories of selected individual electrons (color indicates instantaneous Lorentz factor; the color 
scale is the same as in Figure 5), together with the number density contours (cf. the upper left panel of Figure 2). Arrows indicate the 
main acceleration episodes. The first row shows examples of particles accelerated by reconnection electric field at magnetic X-points. The 
second row shows examples of particles accelerated by anti-reconnection electric field between merging plasmoids. The third row shows 
examples of particles accelerated within the plasmoids. The fourth row shows other cases, including combinations of different mechanisms, 
and possibly acceleration at the post-merger disturbances. 


that seen for particle 128 (representing acceleration at 
the magnetic X-points), which bounces off a plasmoid 
without additional acceleration. We observe additional 
structures in the electric field propagating outward af¬ 
ter each plasmoid merger; these appear to be triggered 
by the merger, carrying excess momentum of the merg¬ 
ing plasmoids, possibly via their Poynting flux (see the 
lower left panel of Figure 2). We associate the sharp ac¬ 
celeration of particle 113 with its interaction with one of 
those disturbances. 


3.2.4. Statistics of the main acceleration sites 

In order to evaluate the relative importance of the ac¬ 
celeration mechanisms discussed in the previous subsec¬ 
tion, we have compared the parameters of all main ac¬ 
celeration episodes for our representative sample of en¬ 
ergetic particles. The lower left panel of Figure 2 shows 
particle locations in the spacetime diagram correspond¬ 
ing to the beginning of each main acceleration episode. 
The color scale indicates the mean value of the z velocity 
(v z ) /c during the main acceleration episode, which is a 
measure of the dominant orientation of the electric field 
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Figure 5. Full 6-dimensional phase space and electromagnetic fields along the trajectories of selected individual electrons. The color 
scale used on the first 4 panels (from the left) and shown on the right indicates the electron Lorentz factor. On the last 3 panels, the 
x components are orange , y components green , and z components magenta. Arrows indicate the main acceleration episodes. Horizontal 
dashed lines indicate specific simulation times discussed for each particle in the main text. 
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accelerating the particle. The spacetime distribution of 
particles undergoing significant acceleration is very com¬ 
plex. We observe several regions in the spacetime dia¬ 
gram where a large number of particles are accelerated. 
These regions are roughly and somewhat subjectively 
classified into three main types defined on the lower right 
panel of Figure 2 : ( 1 ) particle diffusion regions around 
the magnetic X-points where (v z ) > 0.5 c (blue), ( 2 ) re¬ 
gions between merging plasmoids where (v z ) < —0.5 c 
(red), and (3) trailing edges of fast-propagating plas¬ 
moids where (v z ) > 0 (green). Less important accelera¬ 
tion sites are associated with small secondary plasmoids 
in the later stages of the simulation (uot > 1200 ; blue), 
and also merger-triggered disturbances in the early stages 
of the simulation (u; 0 t < 1200). In general, the acceler¬ 
ation of particles begins at uj^t ~ 200 , and is complete 
by ujot ~ 2000 . And while it is widely spread in the 
spacetime diagram, it is not uniform. Very few particles 
are accelerated in the vicinity of static plasmoids, while 
the fast-propagating plasmoids accelerate many particles 
both at the leading and trailing edges (note the avoid¬ 
ance of dense plasmoid cores). 

In Figure 6 , we present the distribution of maximum 
particle energies for the representative sample of ener¬ 
getic particles in the energy range 10 < 7 max < 150. We 
also show separately the contributions from each class of 
the main acceleration sites, calculated for 2 values of pa¬ 
rameter k. For each particle in the sample, we consider 
all m main acceleration episodes, with each episode clas¬ 
sified according to the initial location of the particle in 
the spacetime diagram (Figure 2 , bottom right panel). 
Each episode contributes with weight 1/m to the energy 
histogram at the maximum particle energy 7 max . For 
a consistency check, we compare the histogram of 7max 
values for the entire sample of energetic particles (gray 
points in Figure 6 ) with the momentum distribution of 
all background electrons at the end of simulation (the 
magenta line). We find good agreement between the two 
distributions. The effect of changing parameter k be¬ 
tween 0.5 and 0.25 is modest and limited to 7 max < 40. 

For each class of the main acceleration sites, we per¬ 
formed a fit to the energy distribution of energetic elec¬ 
trons for 7 > 10 with the following power-law model 
N( 7 ) = 7Vo7 _a exp[—( 7 / 70 )^], with either exponential 
(q = 1) or super-exponential (q = 2) cut-off. This ap¬ 
proach is simpler than modeling using both types of cut¬ 
offs adopted in Werner et al. (2014), as our simulation is 
in the region of the (L, oq) parameter space where nei¬ 
ther cut-off shape is expected to dominate (see Section 
4). The best-fit parameters for k = 0.5 and k = 0.25 
are listed in Table 1, and all best-fit functions are shown 
in Figure 6 . In the cases of mergers and other acceler¬ 
ation sites, the exponential cut-off provides a better fit, 
and in the cases of magnetic X-points and plasmoids the 
super-exponential cut-off is preferred. The values of the 
cut-off Lorentz factor are y c ^ 9 — 15 for q = 1 and 
7 C ^ 29 — 36 for g = 2, regardless of the value of k. Par¬ 
ticles accelerated in magnetic X-points have the hardest 
distribution (lowest a), followed by particles accelerated 
in mergers, within plasmoids and at other sites. The 
value of a increases by 0.7 — 1.2 between the fits with 
q = 1 and q = 2 ; since both types of cut-off can be ac¬ 
ceptable, we cannot make a definite measurement of a 
power-law index of the electron energy distribution (cf. 


X-point 1—•—1 plasmoids ■—•—* total 

mergers ■—•—■ other ■—•—■ 



Tmax 


X-point ■—•—« plasmoids ■—•—■ total 

mergers 1 —•— 1 other ■—•—■ 



Figure 6. Distribution of maximum particle energies for the sam¬ 
ple of tracked energetic electrons (black), divided according to our 
classification of main acceleration sites with k = 0.5 (top panel) 
and k = 0.25 (bottom panel): magnetic X-points (blue), regions be¬ 
tween merging plasmoids (red), plasmoids (green), and other sites 
(gray). The vertical scale indicates the actual number of tracked 
electrons in each energy bin. Error bars are Poissonian only. Best- 
fit spectral models - power-laws with exponential (q = 1) and 
super-exponential (q = 2 ) cut-offs (the high-energy tails of the 
q = 1 fits are always to the right from the corresponding q = 2 fits) 
- are shown with dashed lines (see Table 1). The solid magenta 
line shows a normalized momentum distribution of all background 
electrons at the end of the simulation. 

Werner et al. 2014). We find that in the energy range 
10 < 7 max < 20 , the largest contribution (^ 26%) to 
the total energy distribution is from the plasmoids, and 
in the energy range 40 < 7 max < 100, mergers clearly 
dominate (up to ^ 50%) over other acceleration sites. 
The contribution of magnetic X-points is no more than 
~ 25%, and the unclassified (‘other’) acceleration sites 
contribute between ~ 10% at 7 max ^ 70 to ^ 40% at 

7max 00 10- 

In Table 2, we report the statistics of the number and 
classification of the main acceleration episodes for the 
representative sample of 14169 energetic particles with 
7max > 10, and for a representative subsample of 1817 
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Table 1 

Parameters of function IV (7) = No^y~ a exp[— (7/7c) g ] fitted in the range 10 < 7 < 150 to the energy distributions of energetic particles 
according to their acceleration site (see Figure 6). In order to calculate x 2 ; we a< lcl unity to the fitted function; the number of degrees of 

freedom is 12. 




exponential (q = 1) 


super-exponential (q = 

2) 

k 


a 

7c 

x 2 

a 

7c 

x 2 


X-points 

-0.61 ±0.25 

8.9 ± 1.0 

14 

0.61 ±0.15 

29.3 ±2.0 

31 


mergers 

-0.05 ±0.12 

11.9 ±0.8 

10 

0.89 ±0.08 

34.9 ± 1.5 

53 

0.5 

plasmoids 

0.40 ±0.18 

10.9 ± 1.2 

15 

1.31 ±0.08 

31.1 ± 1.4 

11 


other 

1.57 ±0.15 

14.9 ±2.0 

16 

2.21 ±0.08 

35.7 ±2.6 

11 


total 

0.64 ±0.12 

13.2 ± 1.1 

33 

1.43 ±0.06 

35.2 ± 1.4 

85 


X-points 

-0.66 ±0.23 

9.0 ±0.9 

10 

0.55 ±0.13 

29.4 ± 1.7 

22 


mergers 

0.06 ±0.11 

12.8 ±0.8 

10 

0.94 ±0.06 

36.4 ± 1.3 

24 

0.25 

plasmoids 

0.45 ±0.15 

11.9 ± 1.1 

17 

1.31 ±0.06 

33.2 ± 1.4 

12 


other 

1.07 ±0.14 

12.3 ± 1.2 

12 

1.85 ±0.07 

32.5 ± 1.9 

21 


total 

0.65 ±0.12 

13.2 ± 1.1 

33 

1.43 ±0.06 

35.3 ± 1.4 

84 


Table 2 

Statistics of the number and classification of main acceleration 
episodes for different values of parameter k and particle energy 
treshold 7 m i n - For particles with 2 main acceleration episodes we 
report the most frequent sequences of the acceleration sites. Only 
values exceeding 3% are reported. 


k 

0.5 


0.33 

0.25 

7max 

> 10 

> 30 

> 30 

> 30 

1 episode 

38% 

57% 

< 3% 

0% 

X-point 

8% 

14% 

— 

— 

merger 

14% 

29% 

— 

— 

plasmoid 

3% 

5% 

— 

— 

other 

13% 

10% 

— 

— 

2 episodes 

38% 

38% 

50% 

7 % 

X-point —> merger 

— 

3% 

4% 

— 

X-point —> plasmoid 

7% 

12% 

— 

— 

X-point —> other 

4% 

5% 

7% 

— 

merger —> merger 

3% 

4% 

14% 

4% 

merger —>• plasmoid 

— 

4% 

4% 

— 

plasmoid —>• plasmoid 

— 

3% 

— 

— 

other —> merger 

3% 

— 

4% 

— 

other —> plasmoid 

6% 

4% 

— 

— 

other other 

4% 

— 

4% 


3 episodes 

15% 

5% 

33% 

43% 

> 3 episodes 

8% r 

-0% 

14% 

49% 


particles with 7max > 30 for three values of parameter 
k = 0.5, 0.33, 0.25. In the case of k — 0.5, we find that 
the majority of energetic particles with 7 max > 10 were 
accelerated either once or twice (38% each). The most 
important acceleration scenarios are: (1) single acceler¬ 
ation between merging plasmoids (14%); (2) single ac¬ 
celeration in a magnetic X-point (14%); (3) acceleration 
in a magnetic X-point followed by acceleration in a plas- 
moid (9%). The more energetic particles with 7 max > 30 
are more likely to be accelerated only once (57%), partic¬ 
ularly in a single merger event (29%). This is consistent 
with the relative contribution of mergers to the total en¬ 
ergy distribution increasing with the particle maximum 
energy. Two conclusions can be drawn from this: (1) 
while magnetic X-points and regions between merging 
plasmoids are the dominant acceleration sites, relatively 
few particles are accelerated in both of them; and (2) the 
particles that are accelerated within the plasmoids most 
likely were previously accelerated in a magnetic X-point. 


This may help to explain why the contributions to the 
high-energy spectrum from magnetic X-points and plas¬ 
moids are comparable to each other, and different from 
the contribution from mergers. For smaller values of &, 
almost all particles are characterized by multiple acceler¬ 
ation episodes. The most frequent sequence of two accel¬ 
eration episodes is a sequence of two plasmoid mergers. 

4. DISCUSSION 

Relativistic magnetic reconnection is a very complex 
phenomenon, where dissipation of magnetic energy, par¬ 
ticle acceleration, and, under some conditions, high- 
energy emission proceed simultaneously. Investigating 
in detail the distribution and relative importance of the 
main acceleration sites is of the greatest importance for 
understanding the distribution of high-energy radiation 
emitted by these particles. This radiation may be re¬ 
sponsible for a variety of high-energy astrophysical phe¬ 
nomena including gamma-ray flares of blazars and the 
Crab Nebula, gamma-ray bursts, etc. The characteris¬ 
tics of the radiation from plasmoid-dominated relativis¬ 
tic magnetic reconnection will be discussed in detail in a 
separate follow-up work. 

By analyzing a representative sample of energetic par¬ 
ticles, we have provided a comprehensive picture of par¬ 
ticle acceleration for one simulation of relativistic recon¬ 
nection. Other studies that have focused on investigat¬ 
ing the histories of ‘representative’ individual particles 
inevitably introduce subjective choices, which we have 
endeavored to avoid. While we do show representative 
examples for illustrative purposes (Figures 4 and 5), we 
recognize that most individual particle histories are dif¬ 
ficult to classify straightforwardly. For our main results, 
we rely on a simple concept of the main acceleration 
episodes, which can be easily identified for every parti¬ 
cle. Other classification schemes could certainly be con¬ 
sidered. 

Our results indicate that in the plasmoid-dominated 
reconnection regime, there are three classes of main ac¬ 
celeration sites, together accounting for ^ 80% of the 
main acceleration episodes: magnetic X-points, regions 
between merging plasmoids (leading edges of fast-moving 
plasmoids), and trailing edges of fast-moving plasmoids. 
That particles can be efficiently accelerated in relativistic 
magnetic X-points has been demonstrated previously in 
smaller simulations that did not capture multiple plas- 
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moids (Zenitani & Hoshino 2001; Lyubarsky & Liverts 
2008), as well as in previous large simulations (e.g., 
Jaroschek et al. 2004; Zenitani & Hoshino 2005a, 2007). 
That particles can be accelerated between merging plas- 
moids was demonstrated in the non-relativistic case (Oka 
et al. 2010), and mentioned briefly in the relativistic case 
by Sironi & Spitkovsky (2014). Bessho & Bhattacharjee 
(2012) argued that in addition to such acceleration, a 
significant number of particles are decelerated there. We 
have not looked for evidence of such deceleration, how¬ 
ever, we do not think that it would be in contradiction 
with our results. As for the particles accelerated within 
plasmoids, they gain energy mostly when the plasmoids 
attain significant bulk motions, and do so mostly at the 
trailing edge. This acceleration mechanism appears con¬ 
sistent with curvature drift (e.g., Dahlin et al. 2014; Li 
et al. 2015), the importance of which was recently em¬ 
phasized, in the context of relativistic reconnection, by 
Guo et al. (2014, 2015). 

By finding all main acceleration episodes for a repre¬ 
sentative sample of energetic particles, we are able to 
identify the most important acceleration scenarios. Ac¬ 
cording to Table 2, in the case of k = 0.5, the two most 
important acceleration sites are magnetic X-points and 
plasmoid mergers. In the case of k < 0.5, most par¬ 
ticles are accelerated in multiple acceleration episodes, 
and the most frequent sequence of acceleration sites is 
a sequence of 2 mergers. We find that only ^ 4% of 
particles are characterized by a sequence of two main ac¬ 
celeration sites, the first an X-point, and the second a 
merger. We have verified that while the fraction of par¬ 
ticles characterized by this particular sequence increases 
with decreasing value of parameter k in our definition 
of the main acceleration episode, in any case there are 
several other sequences that are more important than 
this one. On the other hand, Sironi & Spitkovsky (2014) 
found - based on a much larger simulation - that most 
of the highest-energy particles are energized in such se¬ 
quences. We suggest that the fraction of particles accel¬ 
erated in plasmoid mergers that previously pass through 
a magnetic X-point may depend on the simulation size. 

We do not find evidence of significant particle acceler¬ 
ation by a Fermi-type process operating in contracting 
plasmoids, as proposed by Drake et al. (2006) in the case 
of non-relativistic reconnection in electron-ion plasma. 
In that mechanism, an electron can be energized when 
circulating around the center of an elongated plasmoid 
by electric fields induced by a gradual contraction of the 
plasmoid. In our simulation, the plasmoids formed in 
mergers are observed to relax on the dynamical time 
scale (i.e., Alfven wave crossing time scale) and subse¬ 
quently they do not undergo significant contraction. On 
the other hand, the first generation of plasmoids form¬ 
ing directly from the tearing mode instability is charac¬ 
terized by significant elongation, as can be seen in Fig¬ 
ure 1 for ujot = 249. These plasmoids circularize before 
ujot = 497; however, few particles are significantly accel¬ 
erated before this time (see the brown line in Figure 3), 
and most of those are accelerated in magnetic X-points. 
We observe many energetic particles circulating within 
plasmoids — a good example is ‘particle 95’ shown in 
Figure 5. However, these particles are significantly ac¬ 
celerated only when the plasmoids they are trapped in 
are themselves accelerated in bulk, as indicated by the 


spacetime distribution of the main acceleration episodes 
in the bottom right panel of Figure 2. We note that the 
evolution of plasmoids proceeds much faster in the rela¬ 
tivistic case than in the non-relativistic case (since in the 
relativistic case, the Alfven velocity is of the order of the 
speed of light). Moreover, the initial episode of contract¬ 
ing islands is only a transient phenomenon, and hence 
this acceleration mechanism is less important in the case 
of relativistic reconnection relative to other mechanisms 
that depend on bulk plasmoid acceleration and mergers. 
Plasma composition may also be important to this prob¬ 
lem, as in pair plasma there is no separation of energy 
scales, and hence the Alfven velocities, between different 
particle species. 

Recent detailed studies of the total particle energy dis¬ 
tribution resulting from relativistic reconnection are re¬ 
ported in Sironi & Spitkovsky (2014); Guo et al. (2014); 
Werner et al. (2014) and Guo et al. (2015). The analysis 
of Werner et al. (2014) suggests that the general form of 
the distribution is a power-law with an exponential [oc 
exp(— 7/7 c )] and super-exponential [oc exp(— ( 7 / 7 c ) 2 )] 
cut-off, depending on the simulation domain size (an ex¬ 
ponential cut-off is expected for simulation domains with 
L x (y)/ cr oPo > 40). The power-law index a (dTV/dy oc 
7 -a ) depends on the upstream magnetization cfq (see 
also Sironi & Spitkovsky 2014; Guo et al. 2014, 2015). 
For > 10, the distribution becomes very hard (a < 2), 
with most energy contained around the high-energy cut¬ 
off which scales linearly with a or L. In the case stud¬ 
ied here — L x = L y = 800po and ao = 16, hence 
L x ( y )/cropo = 50 — we expect to be at the boundary 
between the two regimes characterized by the predom¬ 
inance of the exponential or super-exponential cut-offs. 
Indeed, we find that both types of cut-off can describe 
contributions to the energy distributions of energetic par¬ 
ticles from different types of main acceleration sites (see 
Figure 6 and Table 1). We also find that the compensated 
momentum distribution u 2 N(u) of all background elec¬ 
trons peaks between 10 < ^ pea k < 20 at late stages of the 
simulation, with the effective power-law index a ^ 1.6 
(see Figure 3). 

In principle, the energy distribution of particles could 
extend to even higher energies in a steeper power-law 
tail if at least some particles could undergo multiple 
energizations reminiscent of a first-order Fermi process 
(Hoshino 2012). Our analysis reveals that about 10% of 
energetic particles experience more than 2 main accel¬ 
eration episodes. However, due to our definition of an 
acceleration episode (Ay ^ &7max with k = 0.5), those 
particles must undergo significant deceleration between 
every two accelerations, and the history of those parti¬ 
cles is not consistent with the canonical Fermi process 

(At ~ i)- 

Previous simulations performed - unlike here - in the 
radiatively efficient regime, where the synchrotron cool¬ 
ing time scale for the most energetic particles is compa¬ 
rable to or shorter than the dynamical time scale, em¬ 
phasized particle acceleration in magnetic X-points on 
Speiser-type orbits as the most important acceleration 
mechanism (Cerutti et al. 2013, 2014). In the radia¬ 
tively efficient regime, particle acceleration is limited in 
the presence of a strong perpendicular magnetic field. 
These works demonstrated that particles accelerated in 
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magnetic X-points experience only weak perpendicular 
magnetic field components. On the other hand, one 
can expect stronger perpendicular magnetic field com¬ 
ponents in the plasmoids. Hence, in the radiatively effi¬ 
cient regime acceleration within plasmoids or in plasmoid 
mergers may be suppressed by the severe radiative losses. 
This issue requires a further, more detailed investigation. 

Finally, the simulation presented in this work uses no 
initial magnetic perturbation and periodic boundary con¬ 
ditions, which leads to a rapid hierarchical evolution of 
the plasmoids via balanced mergers that tend to can¬ 
cel the bulk of the x-component of momentum. The 
distribution of the particle acceleration sites may be 
fundamentally different with open boundary conditions. 
Such simulations tend to develop a roughly statistically 
steady-state configuration (with some tearing) with a 
central magnetic X-point and broad Alfvenic outflows 
(Daughton & Karimabadi 2007; Loureiro et al. 2012). We 
would then expect a significantly smaller role of plasmoid 
mergers in accelerating particles and producing sharp ra¬ 
diation flares. 

5. CONCLUSIONS 

We presented a study of the evolution of a plasmoid- 
dominated reconnection layer and of the distribution and 
relative importance of the main particle acceleration sites 
by analyzing a kinetic numerical simulation of relativistic 
magnetic reconnection in pair plasma with initial magne¬ 
tization a o = 16 and no guide field. Substantial insight is 
gained by analyzing spacetime diagrams for various phys¬ 
ical parameters, and by investigating a representative 
sample of tracked energetic particles (see Figure 2). We 
find that the plasmoids develop a robust intrinsic struc¬ 
ture with colder dense cores and hotter outer shells that 
is preserved through multiple plasmoid mergers. After 
the initial episode of elongated plasmoids arising from the 
tearing mode instability, the plasmoids are not contract¬ 
ing noticeably. We identify three major types of particle 
acceleration sites: (1) magnetic X-points with the recon¬ 
nection electric field E z < 0 (the sign refers to this partic¬ 
ular case), (2) regions between merging plasmoids with 
the anti-reconnection electric field E z > 0, and (3) the 
trailing edges of accelerating plasmoids. Notably, parti¬ 
cle acceleration is not observed around static plasmoids, 
and it cannot be associated with the initial transient 
phase of plasmoid contraction. We introduce the concept 
of the main acceleration episode as the shortest period of 
time [ti : t 2 ] where j(t 2 ) - 7 O > Mlmax ~ 7min), and 
use it to compare the relative importance of the different 
classes of acceleration site in producing the high-energy 
electron distribution. For k = 0.5, we find that plasmoid 
mergers are able to accelerate particles to higher energies 
(7 ~ 100), with a harder distribution, than magnetic X- 
points or accelerating plasmoids (the power-law index 
differs by ~ 0.35, but its actual value depends on the 
assumed shape of the high-energy cut-off; see Figure 6 
and Table 1). As a consequence, particles accelerated in 
plasmoid mergers dominate the energy distribution for 
7 > 40. The particles accelerated during plasmoid merg¬ 
ers are usually not the same particles that are accelerated 
at other sites (see Table 2). The particles accelerated in 
accelerating plasmoids were usually also accelerated in 
magnetic X-points. These results will be essential for 
interpreting the time-, angle- and energy-dependent ra¬ 


diative signatures of relativistic magnetic reconnection, 
which will be described in a future study. 
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